Roles of axial anomaly on neutral quark matter with color superconducting phase 
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We investigate eflects of the axial anomaly term with a chiral-diquark coupling on the phase di- 
agram within a two-plus-one-flavor Nambu-Jona-Lasinio (NJL) model under the charge-neutrality 
and /^-equilibrium constraints. We find that when such constraints are imposed, the new anomaly 
term plays a quite similar role as the vector interaction does on the phase diagram, which the present 
authors clarified in a previous work. Thus, there appear several types of phase structures with mul- 
tiple critical points at low temperature T, although the phase diagrams with intermediate-T critical 
point(s) are never realized without these constraints even within the same model Lagrangian. This 
' drastic change is attributed to an enhanced interplay between the chiral and diquark condensates 

(N ■ due to the anomaly term at finite temperature; the u-d diquark coupling is strengthened by the 

relatively large chiral condensate of the strange quark through the anomaly term, which in turn 
' definitely leads to the abnormal behavior of the diquark condensate at finite T, inherent to the 

asymmetric quark matter. We note that the critical point from which the crossover region extends 
to zero temperature appears only when the strength of the vector interaction is larger than a critical 
■ value. We also show that the chromomagnetic instability of the neutral asymmetric homogenous 

two-fiavor color superconducting(2CSC) phase is suppressed and can be even completely cured by 
the enhanced diquark coupling due to the anomaly term and/or by the vector interaction. 
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I. INTRODUCTION 

> . 

: It is generally believed that the strongly interacting matter exhibits a rich phase structure in extreme environment 
\^ such as at high temperature and high baryon chemical potential. Experimentally, RHIC (Relativistic Heavy Ion 
Collider) and LHC (Large Hadron Collider) may provide more information on this topics. Theoretically, some results 
CO have been already obtained on a sound basis: First, the lattice simulations of quantum chromodynamics (QCD) 
\ indicate that, for physical quark masses, the transition from the hadronic phase to the quark gluon plasma (QGP) is a 
. smooth crossover at finite temperature and vanishing baryon chemical potential [H Q , whereas in the low temperature 
and very high density region, the techniques of perturbation QCD can be used and the color flavor locking (CFL) Q 
phase is proved to be the ground state of QCD [4,-^]- 

However, the above methods based on the first principle fail at the low temperature and moderate density region, 
due to the sign problem or the non-perturbative effect. Phenomenologically, such a region in the T-/i plane is more 
relevant to reality and hence interesting since it is directly related to the physics of compact stars. On that account, 
chiral models of QCD such as the NJL model (sl-flU that embody the basic low-energy characteristics of QCD such as 
symmetry properties have been extensively used to explore the T-/i phase diagram of strongly interacting matter. In 
particular, such model calculations suggest that CSC phase may occur at low temperature and large chemical potential 
(for reviews, see [l^ - |T5| ). In addition, a popular result from the model studies is that the chiral phase transition 
always keeps first order at the low-temperature region [l6l - l2l| . Combined with the crossover transition confirmed 
by lattice QCD at zero baryon chemical potential, usually, a schematic T-^l phase diagram with one chiral critical 
point(CP) is widely adopted in the literature [13]. Such a CP may be located at relatively high temperature and 
low baryon chemical potential, which has attracted considerable attention as it is potentially detectable in heavy-ion 
experiments [Hill]. 

Generally, there is no reason to rule out the possibility that the QCD phase diagram may contain more than one 
chiral CP, especially when the chiral and diquark condensates are considered simultaneously; some rich structures 
with multiple CP's may be expected for the phase diagram owing to somehow enhanced interplay between the two 
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types of condensates^. It has been shown on the basis of the NJL model that it is indeed the case [27h3(J |: the QCD 
phase diagram can admit muhiple CP's when the repulsive vector interaction [23, HI] or the charge-neutrality and 
/3-equilibrium [2^ or both of these two ingredients [SOI are included ^ . This is because the two ingredients act so as to 
enhance the competition between the chiral and diquark condensates and thus the would-be first-order boundary line 
in the low-temperature region of the T-fj, plane can be turned into a smooth crossover or multiply-cut crossover lines 
with new CP('s). Indeed it has been found [13] that the number of the chiral CP's may vary from zero to four with 
the joint effect of these two ingredients. Moreover, the present authors have shown [30] that the vector interaction 
can effectively suppress the chromomagnetic instability [33| in the asymmetric homogeneous CSC phase. 

It is noteworthy that a direct coupling term between the chiral and diquark condensates can be supplied by the 
axial anomaly [^, [U, [s^, which thus might lead to a new CP in the low-temperature region, as first conjectured 
in [3^ on the basis of an analysis using the Ginzburg-Landau (GL) theory in the chiral limit; see also subsequent 
detailed analyses [13, US'] though still in the chiral limit. It is, however, to be noted that the GL theory assumes that 
both the diquark and chiral condensates are sufficiently small around the phase boundaries provided that the phase 
transitions are of second-order. In addition, the GL theory itself can not determine the coefficients in the action, and 
some microscopic model or theory is necessary for such a determination. 

Recently, a microscopic calculation has been done [39] with use of a three-flavor NJL model for massive quarks 
incorporating the axial anomaly term with a form of a six-quark interaction [3, [s^ . [ssl [37| , the coupling constant of 
which is denoted by K': It was claimed in [39] that the low-temperature CP can exist owing to the axial anomaly 
for an appropriate range of the model parameters even off the chiral limit but still with a flavor symmetry as in 
the GL approach in the chiral limit. It should be noticed, however, that the SU(3)-flavor symmetry may lead to a 
special type of CSC phase, i.e., the CFL phase, as is taken for granted in [s^ - lsgj . which automatically satisfies the 
charge-neutrality and /3-equilibrium constraints. 

Then one may suspect that the possible emergence of the new CP might be an artifact of such an ideal situation with 
the three-flavor symmetry. Nevertheless, it is a very interesting possibility that the axial anomaly-induced interplay 
between the chiral and diquark condensates would lead to a new CP in the low-temperature region. Thus it is worth 
exploring to see whether a new low-temperature CP is induced by the axial anomaly in a dynamical model of QCD 
by considering the realistic situation with the broken flavor symmetry by the hierarchical current quark masses. We 
note that once the quark mass difference of different flavors is taken into account, it becomes a complicated dynamical 
problem to make the charge-neutrality and /3-equilibrium constraints satisfied. 

More recently, such a realistic calculation in the framework of a two-plus-one-flavor NJL model has been done by 
Easier et al. [40|; they have shown that such a new low-temperature CP is not found in such a model even with 
the axial anomaly term, because an unusual interplay between the chiral and diquark condensates induced by the 
anomaly term actually leads favorably to the 2CSC phase \4U \4 3 rather than the CFL phase near the chiral phase 
boundary, even in the case with the equal quark mass limit |40[ . 

It is worth emphasizin g h ere that the constraints by the charge-neutrality and /3-equilibrium are not taken into 
consideration in Q nor |40j in contrast to [l^, [13] where various types of multiple-CP structures are found in the 
phase diagram. Thus the following two questions arise naturally: Will the results found in [ioj be altered or not when 
the charge- neutrality and /3-equilibrium constraints and/or the vector interaction are taken into account? Or will the 
phase structure with multiple CP's found in 30] rather persist when taking into account the coupling between the 
chiral and diquark condensates induced by such a six-quark interaction? The main purpose of this paper is to answer 
these questions by incorporating the anomaly term that breaks the C/4(l) symmetry as well as the vector interaction 
under the constraints of the charge-neutrality and /3-equilibrium in the two-plus-one-flavor NJL model. The present 
work may be regarded as either an extension of the paper [i^] by incorporating the charge-neutrality, /3-equilibrium 
and the vector interaction, or an extension of the paper [30j by including the ^fC'-interaction. 

The main conclusion we reach is that the key results on the phase structure obtained in Ref. [29ll30j persist even when 
the attractive i^T'-interaction is incorporated. That is, there appear new CP('s) at the intermediate temperature owing 
to charge-neutrality constraint and then the transition in the low temperature region extending to zero T becomes a 
crossover when the strength of the vector interaction becomes larger than a critical value: Thus the number of the 
CP's can be even more than two, depending on the values of some related coupling constants. Strikingly enough, 
we find that the interplay between the chiral and the diquark condensates induced by the anomaly term even acts 
toward realizing the multi-CP structure of the phase diagram under the neutrality and /3-equilibrium constraints even 
without the help of the vector interaction. Accordingly, the results in Ref. j4(|] are modified by considering these 
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constraints. Even though the chiral boundary in the low-T region extending to zero T also remains first order in our 
case in the absence of the vector interaction, which shrinks and vanishes eventually as K' becomes large and exceeds 
a critical value. 

We shall also examine the chromomagnetic (in) stability under the influence of the axial anomaly, as was done in 
[SOi] . It is well known that the asymmetric homogenous 2CSC phase suffers from the chromomagnetic instability. 
At zero temperature, the calculation based on the hard-dense-loop (HDL) method 33] suggests that the Meissner 
mass squared of the 8th gluon becomes negative for ^ > 1 while the 4th- 7th gluons acquire negative Meissner 
masses squared for ^ > 1/a/2; here 5)1 and A denote the difference of the chemical potentials of u and d quarks 
and the gap, respectively. Note that i5/Lt is just equal to a half of the electron chemical potential when the vector 
interaction is absent, and this quantity is to be replaced by an effective chemical potential 5jl (see below) when the 
vector interaction is present, as shown in [30| . The instability of the asymmetric homogenous CSC phase should 
imply the existence of a yet unknown but stable phase in this region of the T-^ plane. Candidates of such a stable 
phase include the Larkin-Ovchinnikov-Fulde-Ferrel (LOFF) phase [isl] and gluonic phase |43|. Besides developing the 
possible new phases, the instability problem may also be totally or partially gotten rid of by some other mechanisms. 
For instance, the instability problem becomes less severe simply at finite temperature because the smeared Fermi 
surface relaxes the mismatch of the Fermi spheres of the asymmetric quark matter [45l-l47| . Furthermore, it is known 
that the lar ger the quark mass and the stronger the diquark coupling, more suppressed the instability even at zero 
temperature[48|. Recently, the present authors[30] have shown that the repulsive vector interaction can also resolve 
the instability problem totally or partially. The stability by the vector interaction is realized due to the following 
two ingredients: (1) the density difference between the u and d quarks reduces the mismatch in the effective chemical 
potentials; (2) the nonzero vector interaction suppresses the formation of high density and hence larger quark masses 
than those obtained without the interaction are realized. We shall show that the new anomaly term play a quite 
similar role as the vector interaction and the interplay between the chiral and diquark condensates induced by the 
anomaly term acts toward suppressing the unstable region of the homogeneous 2CSC phase in the T-/i plane; the 
neutral 2CSC phase can become even free from the chromomagnetic instability if K' is larger than a critical value 
K'^, which can be reduced significantly when the vector interaction is incorporated. 

This paper is organized as follows. In Sec. II, the two-plus-one-flavor NJL model with the extended flavor-mixing six- 
quark interaction is introduced under the constraints of the charge-neutrality and /3-equilibrium. The phase diagram 
of the neutral strongly interacting matter with the influence of the axial anomaly is presented in Sec. III. Sec. IV focuses 
on the role of the axial anomaly on the chromomagnetic (in) stability. The conclusion and outlook are given in Sec.V. 



II. NJL MODEL WITH AXIAL ANOMALY AND VECTOR INTERACTION 



A. The model Lagrangian 



We start from the following two-plus-one-flavor NJL model with the vector interaction [30|, |49j and two types of 
six-quark anomaly terms, 



£ = ^ m )^ + £(4) + + + C% 



where m = diagj-(mu, md, rris) denotes the current-quark mass matrix and 



8 r 



i=0 



3 

Gd [(^^75t{i,^V'c)(^c*75t{i,^^) + (^t{i,^^c)(^ctfi|^) 



= -K det [4, (1 + 75) V] + det (1 - 75) ^] 



T E E [ii'^ltlil ± -i,)^c){^tl{l ± 75)V5c)(Vi.(l ± 75)V', 



(1) 

(2) 
(3) 
(4) 
(5) 



Here the four-fermion interactions are all invariant under the U{i)R x [/(3)L-transformation in the flavor space. In 
our notations, the Gell-Mann matrices in flavor (color) space are \{'^'^^ with i = 1, . . . , 8, and Xq^'^^ = a/2/3 l/(c), and 
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m„,d[MeV] 


m,[MeV] 


GsA^ 




A [McV] 


M„,d [MeV] 


5.5 


140.7 


1.835 


12.36 


602.3 


367.7 


MMeV] 


m^[MeV] 


niK [MeV] 


nir,. [MeV] 


m^[MeV] 


Ms [MeV] 


92.4 


135 


497.7 


957.8 


514.8 


549.5 



TABLE I: Model parametrization of two-plus-one-flavor NJL. 



the antisymmetric one is denoted hy with i — 1,2,3 : 

ip) = A^(=\ 4'^^ = Ap\ t^(^^ = A{(^). (6) 

The scalar interaction in C^"^ is responsible for the dynamical chiral symmetry breaking in the vacuum with the 
formation of the chiral condensate, while the vector interaction can be used to investigate the effect of density-density 
interaction on the chiral phase transition [s^l^. In Eqs.® and (O, 4'c stands for C-ip'^ and C — 17072 is the Dirac 
charge conjugation matrix. We remark that the suffix 3 in denotes the channel for the u— d pairing, for example. 
For lower temperature T and large enough baryon chemical potential /i, >C|j*-' leads to the formation of diquark 
condensate in the color-anti-triplet channel [1, |4l|, H^. Besides the four-fermion interactions, Lagrangian ([Ij also 

contains two types of six-quark interactions, and -C^^j: the former is the traditional Kobayashi-Maskawa-'tHooft 
(KMT) interaction [5l|, [s^ and its effect on the phase diagram in T-fi plane is fully examined [ll|, [H, [H, [53l - [55| . 
whereas the latter could be obtained by a Fierz transformation of the former and induces the coupling between the 
chiral and diquark condensates [1, [13, [13, • We remark that both interactions respect the flavor symmetry 
of SU{3)ii X SU{3)l X U{1) while violating the Ua{1) symmetry as mentioned above. The former is responsible 
for accounting for the abnormally large mass of 77' beyond the Weinberg inequality!!^ (in contrast to other pseudo 
Nambu-Goldstone bosons ) in the effective chiral model and can be identified as an induced quark interaction from 
instantons^52i, i57|]. The introduction of the latter to the Lagrangian expands the study of CSC to the six-fermion 
level m. 



B. The model parameters 

The numerical values of some model parameters are given in Table HI In contrast to jsoj , we only consider the case 
with realistic quark masses. The choice of the model parameters is the same as that in (2^, [sO, (all following 
Ref. [11] ), where Gs, the coupling constant for the scalar meson channel, and K, the coupling constant of the KMT 
term, are fixed by the vacuum physical observables (meson masses and decay constants) . We shall work in the isospin 
symmetric limit in two-flavor space with iriu = rrid = 5.5 MeV, and a sharp three-momentum cut-off A is adopted. 

In contrast to Gs and K, no definite observables in the vacuum are available for determining the coupling constants 
Gv, Gd and K' in such a quark model, although we could read off their values from a Fierz transformation of known 
vertices: the coupling constant K' , for instance, can be related to K through the Fierz transformation of the instanton 
vertex, and K' is found to be identical to K [33 . Since we are mainly interested in the roles of K' and Gy on the chiral 
phase transition and the chromornagnetic instability, both these coupling constants are treated as free parameters in 
the present work. Following (39[ . l40j | . we only consider the attractive interaction between the chiral condensate and 
the diquark condensate. Namely, the coupling K' is kept positive. As for the ratio of Gd/Gs, we adopt the standard 
value from Fierz transformation in this paper. Due to the contribution from the KMT interaction, the ratio Gd/Gs 
from Fierz transformation should be 0.95 rather than 0.75 obtained by only considering the four-quark interaction 
[T3 |. Such a choice of the coupling has also been used in Refs. [1^ and [4o|. In the literature, the diquark-diquark 
interaction near the standard value from Fierz transformation is usually called the intermediate coupling. 



3 We remark that the effects of the vector interaction on the baryon-number susceptibihty and the chiral transition in the two-flavor case 
are examined in [sQl and [1611271]. respectively. 
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C. Thermodynamic potential with the constraints of charge-neutrality and /^-equilibrium 

The grand partition function of the NJL model is given by 

Z = e-"^/^ = / D^D^Pe'I''^"'^^+^'^^\ (7) 



where n is the thermodynamic potential density and fi is the quark chemical potential matrix. In general, the quark 
chemical potential matrix /t takes the form jsoj 

/i = ^ - /ieQ + AisTa + MsTs, (8) 

where fj, is the quark chemical potential (i.e. one third of the baryon chemical potential), fie the chemical potential 
associated with the (negative) electric charge, and /X3 and /ig represent the color chemical potentials corresponding to 
the Cartan subalgebra in the SU(3)-color space. The explicit form of the electric charge matrix is Q = diag(|, — -jj — -j)) 
in flavor space, and the color charge matrices are = diag(i, —^,0) and Tg = diag(i, i, — |) in the color space. The 
chemical potentials for the quarks with respective flavor and color charges are listed below: 



gA^e 2' 3^" ' r^y" ' ^r^^ 2''^' 3*^° ' '5* ' 3'^^ 2^'^ 3^ 

2 2 12 12 

3 A^e — 3M8 , A*M = M + 3/^e — 3/^8 , A^bs = A« + 3/^e — 3/ 



Corresponding to the chiral and diquark interactions in Eq. ([T]), we assume that the following condensates are 
formed in the system, namely, the scalar quark- antiquark condensate 

CT. = (10) 

and the scalar diquark condensate 

S^ = {^cil4n^)■ (11) 

In addition, we remark that the quark-number (or baryon-number) density 

P^ = (12) 

has a finite value for finite fi. Note that the indices 1,2 and 3 in Eqs. ([10]) and represent u, d and s quarks, 
respectively, whereas in Eq. (jlip . the indices 1, 2 and 3 stand for the diquark condensate in d-s, s-u and u-d pairing 
channels, respectively. Here we have assumed the condensates and the density are all hornogeneous; the study of the 
phase structure with inhomogeneous condensates and/or baryon-number density psl I60l463j | is surely intriguing but 
beyond the scope of the present work. 

The constituent quark masses and the dynamical Majarona masses are expressed in terms of these condensates as 
follows: 

K' 

M, = m, - ^Gs(y^ K\e^k\o30k + , (13) 

and 

A, = 2{Gd - ^o,)s, . (14) 
Similarly, it is convenient to define the dynamical quark chemical potential for flavor i by 

fli = fi^- AGvPi ■ (15) 

A few remarks are in order here: 

1. Both types of the anomaly terms -Cj^*-^-' and £^^j contribute to the constituent quark masses in Eq. (fT3|) . and 
thus if K' and the diquark condensate Sj are finite, chiral symmetry is dynamically broken even when the usual 
chiral condensates are absent. 
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2. The new anomaly term also modifies the formula for the Majarona mass for the CSC phase so that the 
chiral condensates affects the Majorana mass, and hence induce an interplay between the two condensates: 
Indeed the 'bare' diquark-diquark coupling Gd is replaced by an effective one, G^. = Gd — ^cn, as shown 
in Eq. (IT4|) . which is dependent on the chiral condensates. Thus the flavor-dependent effective coupling G"^^ is 
now dependent on T and fi through cr^. 

3. Equations. (IT3)) and (IT4|) clearly show that the flavor-mixing occurs not only in the usual chiral condensates due 
to but also in the diquark condensates owing to -C^^j, which would lead to interesting physical consequences. 

4. It is also to be noted that the dynamical quark chemical potential jii for u and d quarks are different from 
each other because of the constraint of electric charge-neutrality (/id > in 2CSC; notice also, however, that 
they are dependent only on the respective density pu,d and hence the dynamical chemical potentials flu^d tend 
to come closer because pd > Pu with the common coupling constant Gy [30l |. 

In the mean field level, the thermodynamic potential for the two-plus-one-flavor NJL with the charge-neutrality 
constraints reads 



n = nl + 2GsY.(T^,-2GvY.pi + Y,{GD~^cr^)\s^ 



2 



2 



- 4ifaia2a3-^5]lndet^. (16) 

p 

Notice the presence of the new cubic-mixing terms among the chiral and diquark condensates. In Eq. (jl6p . fi; denotes 
the contribution from free leptons. Note that VLi should include the contributions from both electrons and muons 
for completeness. Since >> and M,. « 0, ignoring the contribution of muons has little effect on the phase 
structure. Therefore, only electrons are considered in our calculation and the corresponding VLi reads 

Due to the large mass difference between s and u [d] quarks, the most favored phase at low temperature and 
moderate density tends to be the 2CSC rather than CFL phase, as demonstrated in the two-plus-one-flavor NJL 
model [13, [2l|. Surprisingly enough, if the anomaly term C^d'^^^ is incorporated, the 2CSC phase turns to be still 
favored in the intermediate density region even when the three flavors have the equal mass '4(11 . Needless to say, the 
dominance of the 2CSC phase over the CFL one is more robust when the realistic mass hierarchy for the three flavors 
is adopted. Moreover, it is worth mentioning here that the mass disparity favors the 2CSC phase with the u-d pairing 
also through the inequality of the effective diquark coupling G'j-,^ > G'jj^ 2 when the anomaly coupling K' is present; 

see Eq. (HH). Since the main purpose of the present work is to explore how the axial anomaly term C^^^ affect the 
phase boundary involving the chiral transition at moderate densities, under the constraints of the charge-neutrality 
and /3-equilibrium, we only consider the 2CSC phase in the following. Note that p-^ in (jS]) vanishes in the 2CSC phase 
because the color SU(2) symmetry for the red and green quarks are left unbroken. 

The inverse quark-propagator matrix in the Nambu-Gor'kov formalism takes the following form in the mean-field 
approximation, 

with 

[G^]-' ^MiL^n±p)~i-p-M, (19) 

where M — diagj(M„, Md, Ms), p = diagj(/iu, Pd, Ps) and a;„ — {2n + l)nT is the Matsubara frequency. Taking the 
Matsubara sum, the last part of the thermodynamic potential (jl6p is expressed as 

"I^El^det^^-f: / (0{(i^.-i?°)+2Tln(l + e-^^/^)}, (20) 



with the dispersion relations for nine quasi-particles (that is, three flavors x three colors; the spin degeneracy is already 
taken into account in Eq. (1^01) ) and nine quasi- antiparticles. In Eq. (pU)) . represents Ei{M = m, A = 0, p = 0). 
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FIG. 1: The phase diagrams in the T-^ plane for various values of K' in the two-plus-one-flavor NJL model with the charge- 
neutrality and /3-equilibrium being kept. The vector interaction is not taken into account. The thick solid line, thin solid line 
and dashed line denote the first order transition, second order transition and chiral crossover, respectively. 



The s quark and unpaired blue u and d quarks have twelve energy dispersion relations with a similar form. For 
example, the dispersion relations for the blue u quark and anti blue u quark are 

Ebu = E - jlbu and E^u = E + fi^u , (21) 



respectively, with E = + M^. In the rd-gu quark sector with pairing we can find the four dispersion relations, 

(22) 



Erd-nu =Ea± \{i^rd " Mgu) =El^±5jl, 



rd-gu ^ 2 

and the ru-gd sector has another four as 



rd-gu ^ 2 ^ 

^fd-nu ^ Ea± |(Ard - Mffu) ^ EA±Sfj., 



^ru-gd ^ Ea± ^{f^ru " f^gd) = EaT^JI, 



^ru-Qd ^ Ea± |(Aru - Agd) = T , 



(23) 



ru-gd ^ 2 



where Ea = \/ {E — fi)'^ + and Ea — \/ {E -\- /x)^ -I- A^; from now on A stands for A3. The average chemical 
potential is defined by 
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FIG. 2; The temperature dependence of M^, Ms, /S. and 5ji for fixed K' / K = 2.25 and three different chemical potentials 
fj, = 312 MeV, 320 MeV and 330 MeV. The constraints of electric charge-neutrality and /3-equilibrium are imposed while the 
vector interaction is not incorporated. 



and the effective mismatch between the chemical potentials of u and d quarks takes the form 

i(Me-4Gy(p2-pi)). (25) 

Ignoring the the mass difference between u and d quarks, the determinantal term in Eq. (1161) has an analytical 
form which greatly simplifies the numerical calculation. Adopting the variational method, we get the eight non-linear 
coupling equations 

dn _ dVt _ dVt _ dn _ dn _ dn _ dVL _ dD. 
dai 8(73 ds3 dpi dp2 dps dpe ^^ls 

Since is tiny around the chiral transition region pol . [2]| , we shall set it zero with little effect in the numerical 
results [sol. Thus, Eq. ([^^ is then simplified to a system of seven coupled equations. 



III. PHASE STRUCTURE WITH THE AXIAL ANOMALY 



In this section, we show numerical results of the effects of the new six-quark interaction ([5]) on the chiral phase 
transition under the charge-neutrality and /3-equilibrium constraints with or without the vector interaction. Since we 
are mainly interested in the phase diagram involving chiral transition at low temperatures, all the phase diagrams 
will be plotted for the range 250MeV < < 400MeV where the chiral transition is expected to be relevant. As for the 
type of the CSC phase, Ref. ^Qj indicates that the CFL phase is only reaUzed for p > 460MeV even when K' ~ with 
the same model parameters as ours, and an increase of K' pushes the CFL phase to even higher p region. Therefore, 
we exclusively consider the 2CSC phase near the chiral boundary without a loss of generaUty. 

In the following, we use the same notations as in Ref. [29l , ISG . l36j to distinguish the different regions in the T-p 
phase diagram. Namely, NG, CSC, COE, and NOR refer to the hadronic (Nambu-Goldstone) phase with cr ^ and 
A = 0, color-superconducting phase with A ^ and a = 0, coexisting phase with a ^ and A ^ 0, and normal 
phase with ct = A = 0, respectively, though they have exact meanings only in the chiral limit. 



A. The case without vector interaction 



We first show the numerical results in the case without the vector interaction. The phase diagrams with varying 
coupling constant K' are displayed in Fig.[TJ In contrast to Fig. 8 in Ref. [1^, the multi-CP structure can still appear 
with a choice of K' in the phase diagram when the charge-neutrality, /3-equilibrium and the new axial anomaly term 
are simultaneously taken into account. For K' /K = 2.0, Fig.[T^ shows that there exists only one usual chiral CP even 
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though the COE emerges: We remark that the COE region does not exist when K' /K = 0, which is not displayed in 
Fig. [TJ Figure [T|d shows that when K' /K is increased to 2.25, the chiral transition turns to a crossover at relatively 
lower temperatures, and hence there appear two new chiral CP's. With a further increase of K' / K ^ the boundary 
line for first-order transition at higher temperature shrinks and thus the two crossover boundary lines in Fig. [Ud join 
with each other, and eventually only one CP is left in the phase diagram, as shown in Fig. [TJ:. When K' / K is large 
enough. Fig. [TJl indicates that the first order boundary vanishes completely and there is no chiral CP in the phase 
diagram. 

We note that the emergence of the three CP's in Fig. [Ud comes from a joint effect of the interplay between the 
chiral and diquark condensates and the electric charge-neutrality constraint. First of all, we recall that the abnormal 
thermal behavior of the diquark condensate that it has a maximum at a finite temperature in the COE is responsible 
for the emergence of the multiple chiral CP structure [27l - [30} . Such a behavior is also observed in the present case, 
as displayed in Fig. [T]d. As first indicated in Ref. (29j . when /ig = /i^ — /iu is positive, the boundary of the chiral 
transition is shifted towards higher region, and leads to the formation of the COE at low-temperature region, in 
which the chiral phase transition is significantly weakened by the smeariiig of the Fermi surface inherent in the CSC 
phase. In this regard, /ie plays a role of an effective vector interaction f2^, H^l- On the other hand, the chiral 
anomaly term with positive K' intensifies the competition between the chiral and diquark condensates due to the 
enhanced effective diquark-diquark interaction. Thus when K' is increased, the CSC region expands towards lower /i 
region in the T-/i plane. Consequently, the COE region tends to be more easily formed when both /ig and K' take 
effects. Therefore the chiral transition is significantly weakened and the smooth crossover gets to appear with new 
CP's in the intermediate temperature owing to the abnormal thermal behavior of the diquark condensate. 

The T-dependence of M„,Ms, A and 5fi for fixed K' /K = 2.25 and several values of jj, is shown in Fig. [51 One 
can see that, with increasing T, the constituent quark masses decrease persistently while the Majarona mass for CSC 
first increases, has the maximum value and then decreases in the COE region and nearby. Let us see the details for 
each value of /x. For a small ji = 312 MeV, the diquark pairing is weak and the gap A does not appear at lower 
temperature region. Thus the chiral phase transition keeps the nature of the first order at Tci ~ 75 MeV. At a 
lager /i = 320 MeV, the diquark pairing becomes significant and the A shows the abnormal thermal behavior with a 
maximum value around Tc2 ~ 60 MeV, and hence the chiral phase transition turns to a smooth crossover owing to 
the competition with the diquark condensate in the COE region. For even larger fi = 330 MeV, the diquark pairing 
becomes more significant and the A still shows the abnormal thermal behavior. However, the competition between 
the two condensates is not strong enough to qualitatively change the nature of the chiral restoration and a first order 
transition happens at Tcs w 40 MeV. The reason why the crossover does not occur at fj, = 330 MeV but happens at 
fi = 320 MeV can be understood as follows: Starting from the same point (T = Tc3, /i = 320 MeV) in the COE region, 
an increase of T affects the nature of the chiral transition more significantly than that of n does, since Tea < Tc2- 

We have seen that the abnormal thermal behavior of the gap A plays an essential role in realizing the multi-CP 
structure of the phase diagram. Such an unusual T-dependence of the A can be attributed to the following two 
mechanisms: (i) the mismatch between the chemical potentials of u and d quarks owing to the charge-neutrality and 
/3-equilibrium constraints and (ii) the small Fermi spheres of the quarks in the COE region due to the relatively large 
quark masses: First, the difference in the chemical potentials dfl or the mismatch of the Fermi momenta disfavors 
the u-d pairing at zero or small temperature. However, as the temperature is raised in the low-T region, more and 
more u and d quarks tends to participate in the pairing due to the smearing of the Fermi surfaces, especially that of 
the u quark. Of course, when T is raised too much, the pairing will be gradually destroyed. Thus the A will have 
the maximum value at a finite T and then disappears eventually when T is further raised. These dual effects of the 
temperature on the diquark pairing lead to the abnormal behavior of the A. This behavior becomes more prominent 
in the 2CSC phase for a weak diquark coupling [gJ] or in the COE region for a moderate or strong diquark coupling 
[29I l30j . Second, when T is raised, the dynamical quark masses decrease and hence the Fermi spheres or momenta of 
u and d quarks grow significantly for a fixed /it, which means that the density of states at the Fermi surface increases 
with T, and thus the diquark pairing is enhanced in the COE region. The increased diquark condensates in turn tend 
to further suppress the dynamical quark masses. 

Notice that such an increase of the diquark condensate along with increasing T is expected to be most promi- 
nent around the phase boundary of the chiral transition, including the COE region, where the chiral condensates 
change most significantly. For the neutral 2CSC, once the COE is formed, both of these mechanisms take effects 
simultaneously and are mutually enhanced, and thus the formation of the multiple-CP structure is readily made. 

This may explain why no intermediate-temperature CP is realized in Ref. where the chiral-diquark interplay is 
embodied by the anomaly term but without the charge-neutrality and /3-equilibrium constraints; the anomaly term 
solely is insufficient for realizing the abnormal thermal behavior of the A. It should be stressed that the mechanism 
for the emergence of the intermediate-temperature CP's in Fig. [Hd is apparently similar to that in the two-flavor case 
found in [29\. However, the strange quark plays an important role in the present case since the chiral condensate 
of the strange quark contributes positively to the effective diquark-diquark coupling for u and d quarks through 
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FIG. 3: The T-fj, phase diagrams of the two-plus-one-flavor NJL model for several values of K'/K and fixed Gv/Gs = 0.25, 
where the charge-neutrality constraint and /3-equilibrium condition are imposed. With the increase of K'/K, the number of the 
critical points changes and the unstable region characterized by the chromomagnetic instability (bordered by the dash dotted 
line) tends to shrink and ultimately vanishes in the phase diagram. The respective meanings of the various types of lines are 
the same as those in Fig. 1. 



the axial anomaly. We should stress that apart from the appearance of intermediate-temperature CP's, there is a 
common feature with and without charge-neutrality constraint: the chiral transition in the low-T region extending 
zero temperature keeps first order provided that K' does not exceed a critical value at which the first-order line 
completely disappears. 

Last but not least, we remark that the T-^, region where the two new low-temperature CP's are located in Fig.[T]D 
is free from the chromomagnetic instability, which is obvious from Fig. [5^; a detailed discussion on this point will be 
given in Sec. IV. 



B. The case for nonzero vector interaction 



In this subsection, we will investigate the phase diagram when both the vector and the new six-quark interactions 
are present under the charge-neutrality and /3-equilibrium constraints. 

There are some choices for the value of the vector coupling: The chiral instanton-anti-instanton molecule model [H^] 
gives the ratio Gy/Gs = 0.25, while the Fierz transformation of the vertex given in the truncated Dyson-Schwinger 
model [HHl gives the ratio 0.5. Thus we rather treat the Gv/Gs as a free parameter in the range, 0-0.5. 

We first explore the phase diagram in the T-^ plane by varying the ratio K'/K but with Gy/Gs being fixed as 
0.25, the value given in the instanton-anti-instanton molecule model. When K'/K is small and less than 0.5, only 
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the usual phase structure with single CP is obtained. When K' /K exceeds 0.5, other four different types of the CP 
structures appear, as displayed in Fig. [31 At K' /K — 0.55, a phase diagram similar to that in Fig.[T]D is obtained, as 
shown in Fig.[3K, where two new intermediate-temperature CP's emerge. When K' /K is slightly increased to 0.57, the 
chiral transition becomes crossover in the lower-temperature region which extends to zero temperature; thus the total 
number of the CP's becomes four, which indicates stronger competition between the chiral and diquark condensates 
at relatively larger /i. Further increasing K' /K, the low-temperature chiral boundary totally turns into a crossover 
one and only one first-order transition line with two CP's attached remains in the phase diagram, as displayed in 
Fig. Ht. In this case, the number of the CP's is reduced to two accordingly. When K' / K is large enough. Fig. 
shows that only chiral crossover transition exists in the phase diagram with no CP. 

In comparison with Fig. [T] where the vector interaction is not included. Fig. [3] indicates that the phase structures 
with multiple CP's can be realized with relatively small K' owing to the vector interaction. We remark that all the 
types of the chiral CP structures displayed in Fig. [3] by varying K' are obtained by varying Gv without the anomaly 
term js^. In the present case, the number of the critical points changes as 1— >■ 3— >-4— >2— >-0 when K' is increased. 

As is mentioned before, the Fierz transformation of the instanton vertex leads to the identity K' — K, so it is 
of special interest to investigate the phase diagram in the case oi K' = K. A series of phase diagram with fixed 
K' /K = 1 but varied Gy are shown in Fig. S) One finds that all the chiral CP structures in Fig. |3] still appear in the 
phase diagrams, and moreover, even as large as five CP's can exist in the phase diagram, as shown in Fig. [4]:. This 
suggests that the interplay between the chiral and the diquark condensates in the COE region becomes complicated 
once the charge-neutrality, the vector interaction and the axial anomaly are all taken into account. A comparison 
with the case of vanishing K' , which is given in Fig. 7 in [s^l, shows that the parameter range of Gy for realizing 
the low-temperature CP's moves towards lower Gy, which is actually natural because K' gives the same effect as Gy 
on the chiral transition. It is noteworthy that such lower values of Gy are also close to the standard value of Gy / Gs 
derived from the instanton model. The number of the CP's changes asl— 3— >-5— >-4— >-2 — with increasing Gy. 

The anomaly terms in Eq.(IT]) are supposed to originate from the instantons, which are to be screened at finite 
chemical potential and temperature [131 • Accordingly, both the coupling constants K and K' are expected to diminish 
around the phase boundary. However, we emphasize that the main effect of K' is to enhance the chiral condensate 
of the strange quark and the u-d diquark condensate by each other through the cubic coupling among them for the 
realistic quark masses, and Figsl3] and |4] tell us that even smaller values of K' expected at low temperature and 
moderate density can still lead to a quite different phase structure with multiple CP's when the vector interaction is 
present under the charge- neutrality constraint. 



IV. THE INFLUENCE ON THE CHROMOMAGNETIC INSTABILITY 



In this section, we investigate effect of the new axial-anomaly term on the chromomagnetic instability of the 
asymmetric homogeneous 2CSC phase by varying K'. We shall show that the anomaly-induced interplay between 
the chiral and diquark condensates acts toward suppressing the unstable region of the homogeneous 2CSC phase in 
the T-fi plane. Thus the 2CSC phase can become even free from the chromomagnetic instability provided that K' is 
larger than a critical value K'^, which can be reduced significantly when the vector interaction is incorporated. 

The magnetic instability region in the T-fi plane is determined by calculating the Meissner masses squared which 
can be negative when the charge-neutrality constraint is imposed. Here we adopt the same method as that in (45j to 
evaluate the Meissner mass squared 

ml = ^mA)~n{A^O)]B=o, (27) 

where B has the same meaning as that in [45| . Since the strange quark does not take part in the diquark pairing in 
the present case, we can directly use the formula for two-flavor NJL model to calculate the Meissner mass squared. 

The effect of the coupling constant K' on the chromomagnetic instability is shown in FiglSJ We have adopted the 
model parameters in Table U to calculate the Meissner mass squared. Figure [5^ displays the change of the unstable 
region of the chromomagnetic instability with varying K' when the vector interaction is not included. One can see 
that the instability region tends to shrink with increasing K' and eventually vanishes for K' / K > 0.8. This suggests 
that the neutral homogenous 2CSC phase will be totally free from the chromomagnetic instability \i K' = K that 
is derived by the Fierz transformation from the usual instanton vertex. When taking the vector interaction with 
Gy/Gs = 0.5, the unstable region shrinks more significantly with increasing K' and eventually disappears in the 
r-/i plane for K' / K > 0.55, as shown in Figl5}3. This could be an expected result because of the effect of the vector 
interaction on the instability problem found in (soj . 

The reason for the suppression of the chromomagnetic instability by K' is understood as follows. First of all, 
Eq. (|14p tells us that the u-d diquark coupling is enhanced by the presence of the s quark chiral condensate due to 



12 




90 




1 ' 


1 ' 1 ' 


1 ' 


80 






NOR 




70 










60 










^ 50 










|40 




NG 




CSC - 


30 










20 






/ coe\ 




10 

















, / . ,t . 


I 



300 320 340 360 
H[MeV] 

(c) Gv/Gs = 0.195 



380 



400 



90 




' 1 ' 


1 ■ 1 


1 ' 


80 






NOR 




70 
60 










^ 50 
1 40 




NG 




CSC - 


30 










20 










10 






/ coe\ 









1 


1 / . 1 ■ 





300 320 340 360 
|u[MeV] 
(d) Gv/Gs = 0.197 



380 



400 




300 320 340 360 
|j[MeV] 

(e) Gv/Gs = 0.23 



380 



400 




300 320 340 360 
^[MeV] 

(f) Gv/Gs = 0.3 



380 



400 



FIG. 4: The phase diagrams in the two-plus-one-flavor NJL model for fixed K' / K — 1.0 with Gv/Gs being varied, where 
the charge-neutrality constraint and /^-equilibrium condition are taken into account. The meanings of the different line types 
are the same as those in Fig. 1. The number of the critical points changes along with an increase of Gv/Gs- AH the phase 
diagrams are free from the chromomagnetic instability. 
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FIG. 5: The boundary between the stable and unstable homogenous 2CSC regions with (right figure) and without (left figure) 
the vector interaction in two- plus-one-flavor NJL model. With the increase of the ratio K' / K = R, the unstable region with 
the chromomagnetic instability in the T-^ plane shrinks and eventually vanishes. 



the coupling between the chiral and diquark condensates induced by the K' term. On the other hand, as was first 
shown in |48:] through changing the diquark coupling by hand, the chromomagnetic instability tends to be suppressed 
in the strong coupling region and can be completely gotten rid of when the diquark coupling is strong enough. Thus 
we see that the coupling between the u-d diquark and the chiral s-quark condensates by the K' term leads to the 
suppression of the chromomagnetic instability. This is a new mechanism of the stabilization of the gapless 2CSC 
phase, found in the present work. We here emphasize the important role of the strange quark and the anomaly term 
in suppressing the instability: In contrast to the pure two-flavor case, the rather large chiral condensate of the strange 
quark enhances the diquark coupling between the u and d quarks owing to the axial anomaly in the two-plus-one 
flavor case, and this enhancement of the diquark coupling causes the stabilization. 

Due to their common effects on the chromomagnetic (in)stability, Fig[5] suggests that the instability may be totally 
cured in the asymmetric homogeneous 2CSC phase when the coupling constants of the vector interaction and the 
extended six-quark interaction are in an appropriate range. Admittedly, the present work has only dealt with the case 
of the so called intermediate diquark coupling. Nevertheless, for a weaker diquark coupUng, it is expected that the 
system can be still free from the chromomagnetic instability only with larger coupUngs for both the vector interaction 
and the anomaly K'-tevvci. 



V. CONCLUSIONS AND OUTLOOK 



We have explored the phase structure and the chromomagnetic instability of the strongly interacting matter under 
the charge-neutrality constraint within a two-plus-one-flavor NJL model by incorporating a new anomaly term as 
well as the conventional KMT interaction. The anomaly terms have the forms of six-quark interactions and violate 
the Ua{^) symmetry as a reflection of the axial anomaly of QCD. Similarly to the KMT term, the new anomaly 
interaction with the coupling constant K' also induces a flavor-mixing which leads to a direct coupling between the 
chiral and diquark condensates. 

We flrst investigated the role of the axial anomaly on the emergence of the low or intermediate-temperature CP('s) 
without the vector interaction. Owing to the large strange quark mass, the favored CSC phase near the chiral 
boundary is 2CSC rather than CFL, where the electric chemical potential /Zg required by the charge-neutrality plays 
an important role on the chiral phase transition [29| . The once-declared low-temperature CP in the symmetric three- 
flavor limit [1^ was ruled out in Ref. 40] due to the actual dominance of the 2CSC over the CFL. We have shown that 
this is true under charge-neutrality constraint without the vector interaction; the chiral transition in the low-T region 
extending zero temperature keeps flrst order provided that K' does not exceed a critical value at which the flrst-order 
line completely disappears. However, the new chiral anomaly term enhances the competition between the chiral and 
diquark condensates under charge- neutrality constraint, and gives rise to the intermediate-temperature CP's for an 
appropriate range of K' . No such intermediate-temperature CP's had been found in the same model when only either 
the charge-neutrality constraint or the axial anomaly is exclusively included, as shown in [2Q| and ^40]; both of which 
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did not take into account the vector interaction, either. 

We then investigated the T-fj, phase diagram by incor por ating the repulsive vector interaction as weU: We remark 
that this task may be viewed as an extension of the work [30] , in which the effect of the vector interaction on the phase 
diagram is fully explored under the charge-neutrality constraint, to incorporate the anomaly term. We have found 
that the cubic coupling between the chiral and diquark condensates induced by the axial anomaly does not affect 
the qualitative results obtained in jso! ]. Rather, the vector interaction and the anomaly term jointly act so that the 
multiple CP's are realized. Indeed, by varying K' with fixed vector coupling or vise verse, we have shown that all the 
types of multiple-CP structures obtained in 30] can be reproduced. In particular, the phase transition in the low-T 
region extending zero temperature becomes a crossover only when the vector interaction is present with a strength 
larger than a critical value. Furthermore, the maximum number of the CP's can reach as large as five when both 
the interactions are put on. In this case, the low- and intermediate-temperature CP's can appear even with small 
values of K' owing to the help by the vector interaction. This is very welcome because K' in the realistic situation 
at moderate and high density should be weaker than that in the vacuum, since the anomaly term is supposed to 
originate from the instanton configuration which is expected to be suppressed at finite density. 

Besides the influence on the chiral phase transition, we have shown that the axial anomaly also plays an important 
role on the suppression of the chromomagnetic instability for the asymmetric homogenous 2CSC phase, which is 
first disclosed in the present work: With an increase of the extended six-quark interaction, the T-fi region with the 
chromomagnetic instability shrinks and eventually vanishes when the coupling K' is sufficiently large. In particular, 
when taking into account the vector interaction simultaneously, the chromomagnetic instability is suppressed more 
significantly and can be completely gotten rid of by the axial anomaly. 

A general and remarkable message obtained from the present investigation is that the strange quark can significantly 
affect the properties of the neutral strongly interacting matter in which the 2CSC phase with u-d pairing is realized: 
Even though the strange quark does not directly participate in the Cooper pairing in the 2CSC, the interplay between 
the u-d diquark condensate and the strange chiral condensate induced by the anomaly term can lead to drastically 
different phase structure in the T-fi plane under charge-neutrality constraint. 

It should be remarked here that the contribution of other possible cubic flavor-mixing terms composed of different 
condensates, such as 

ap^ ^e'^'^a.p.pk, (28) 
which arise from another type of six-quark interaction 

^e'^"'=e'""(V;,7'^(l±75M)(^j7M(l±75)^r„)(^fe(l±75)V'«), (29) 

are all neglected in Eq. for simplicity. The interaction can be derived from the KMT interaction, which 
may or may not affect the phase structure. Beside their direct contribution to the thermodynamic potential, these 
flavor-mixing terms also modify the dispersion relations of the quasi-quarks: For example, the dynamical quark mass 
becomes dependent on the quark- number density through the term ap^. It is certainly an interesting problem to 
explore the possible effects of these cubic coupling terms on the phase diagram, we leave such a task to a future work. 

Even apart from the neglect of the above vertex (|29p . there are some caveats with the present study based on a 
chiral model that does not embody the confinement effect, and is relied on the mean-field approximation. The results 
obtained in the current study are largely parameter dependent and bears the shortcomings inherent in the mean-field 
approximation. For instance, the result that there can be multiple CP's associated with the chiral transition and 
the CSC actually may merely mean that the QCD matter is very soft for a simultaneous formation of the diquark 
and chiral condensates coupled with the baryonic density along the phase boundary. Of course, a study which 
incorporates these fluctuations should be performed, say, by means of the nonperturbative/functional renormalization 
group method [1^, [G^I with the present model used as a bare model. More profoundly, the effect of the confinement 
should be incorporated even in an effective model approach, which is a more challenging problem since the mechanism 
of confinement is still unclear. Anyway, further studies based on different models and/or methods are needed to 
determine whether the low-temperature CP('s) exists. One of the tasks of future is exploring whether the low- 
temperature CP ('s) persists or not when the inhomogeneous phases are taken into consideration such as the chiral 
crystalline phase [60l463j | or the LOFF phase psj . 
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